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Abstract 

The Brenier optimal map and the Knothe-Rosenblatt rearrangement 
are two instances of a transport map, that is to say a map sending one 
measure onto another. The main interest of the former is that it solves 
the Monge-Kantorovich optimal transport problem, while the latter is very 
easy to compute, being given by an explicit formula. 

A few years ago, Carlier, Galichon, and Santambrogio showed that the 
Knothe rearrangement could be seen as the limit of the Brenier map when 
the quadratic cost degenerates. In this paper, we prove that on the torus 
(to avoid boundary issues), when all the data are smooth, the evolution 
is also smooth, and is entirely determined by a pde for the Kantorovich 
potential (which determines the map), with a subtle initial condition. The 
proof requires the use of the Nash-Moser inverse function theorem. 

This result generalizes the ode discovered by Carlier, Galichon, and 
Santambrogio when one measure is uniform and the other is discrete, and 
could pave to way to new numerical methods for optimal transportation. 
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1 Introduction 



Although optimal transport theory has far-reaching applications, in fields as di- 
verse as continuum mechanics, statistics or image processing, its underlying 
problem is quite simple: how to send one probability measure onto another, while 
minimizing some cost of transportation? Let us denote by fi and v those two mea- 
sures, defined respectively on X and Y. They could for instance represent the 
respective distributions of some goods being produced, and the needs for them, 
and the problem would then be to determine how to organize the supply so that 
the total cost of transportation is as small as possible. 

What we are looking for is a map T : X ^ Y telling us where to send what 
is in x; but T will be suitable only if, for any measurable set A C Y, the goods 
sent by T in A match the needs of the same region, that is to say if fi{T~^ (A)) = 
iy{A). If this condition is satisfied, is said to be the push-forward of fj, by T, and 
we write = T^fj,. Let us denote by c(x, y) the cost for going from x to y, then 
the total cost of transportation we want to minimize is 

/ c{x,T{x))dfi{x). (1) 
Jx 

Notice however that an optimal map may well not exist, and worse, there might 
even be no map transporting /i onto v at all, e.g. if fi is discrete and v is uniform. 

The problem of finding a map T minimizing (1), and such that u = T^fj,, was 
first studied by Monge [1] in the 18th century. In the 1940s, Kantorovich [2] 
introduced the following relaxation of Monge's problem: instead of sending all 
that is in to a unique destination y = T{x), he allowed some splitting. Any 
strategy for sending onto v can then be represented by a measure 7 on X x y, 
such that 7(A x B) gives the share of the goods to be moved from Ax.o B. A 
plan 7 is suitable if it matches the production and the needs, i.e. if 

7(Axy) = MA), ^{XycB) = v{B). 

This simply means that /x and u must be the marginals of 7. Let us denote by 
r(/x, v) the set of all such suitable plans. The total cost of transportation with the 
plan 7 is 

/ c(a;,y)d7(x,y). (2) 
J XxY 

The Monge-Kantorovich problem consists in finding 7 G r(^, v) minimizing (2). 
It is indeed a relaxation of the Monge problem, since if T : X — t- y sends ^ onto 
ly, then the push-forward 7 := (id, T)^fi of ^ by x (x, T(x)) is in r(/i, v), 
and the costs (1) and (2) are equal. 

At the end of the 1980s, Brenier [3, 4] discovered the optimal transport map for 
the Monge problem to exist as the gradient of a convex function and to be unique, 
at least when X = Y = M^, for the cost c{x, y) = jja^ — yp, if is absolutely 
continuous and if and v have finite second order moments. His result was then 
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extended to measures defined on the torus T by Cordero-Erausquin [5], or more 
generally on a Riemannian manifold by McCann [6] . While on the optimal 
map is T(x) = V(/?(x) with (/? convex, on the torus the optimal map can be 
written as T{x) = x — Vil){x) with ip : T'^ — )■ such that ip : x ^ ^x"^ — ip{x) 
defines a convex fiinction on R^. More generally, on a Riemannian manifold 
T{x) = exp^{—V'ilj{x)) for some map ip, called the Kantorovich potential be- 
cause it is linked to a dual formulation for the relaxed problem. 

Being able to compute the optimal map T, or the underlying potential tp, is 
obviously of huge interest. When the measures are discrete, if there is a solution 
to the Monge problem, it can be obtained for instance with the auction algorithm. 
In the continuous case, the solution is also easy to compute in dimension 1, for if 

and I' are absolutely continuous, and if F, G stand for their respective cumu- 
lative distributions, i.e. F{x) := ^((— oo,rE]) and G{y) := cxd, y]), then the 
optimal transport map is T = G^^ o F. 

Unfortunately when the dimension is > 1, there is no such easy formula, and 
it is much more complicated to compute Brenier's map— although not impossible. 
Among the most notable methods, we could cite the one due to Benamou and Bre- 
nier [7], relying on a dynamic formulation of the Monge-Kantorovich problem, 
in which one tries to minimize the average kinetic energy of the particles during 
their transportation. On the other hand, Angenent, Haker, and Tannenbaum [8] 
proposed a steepest descent method, starting from a transport map (for instance, 
the Knothe-Rosenblatt rearrangement) and letting it evolve so as to reduce the 
associated transport cost. A couple of years later, Leoper and Rapetti [9] used 
the characterization of the optimal transport map through the existence of a con- 
vex potential, to compute Brenier's map starting from any potential and, with a 
Newton algorithm, altering it so as to finally get the optimal potential. 

Our hope here is that the results presented in this paper might lead to yet 
another approach for computing Brenier's map. Our starting point is a direct 
connexion, proved by Carlier, Galichon, and Santambrogio [10] a few years ago 
but hinted beforehand by Brenier, between the optimal transport map and the 
Knothe-Rosenblatt rearrangement. This leads us to believe it might be possible 
to compute Brenier's map starting from the rearrangement (as in the paper by 
Angenent, Haker, and Tannenbaum [8]), and then proceeding with a continuation 
method (as in the work by Loeper and Rapetti [9]). 

This so-called Knothe-Rosenblatt rearrangement, which is also built so as to 
send one measure onto another, was first introduced by Rosenblatt [11] and 
Knothe [12]\ It can be defined for absolutely continuous probability measures 
on or on = M^/Z^ (in higher dimension, the construction is analogous) as 
follows: To begin with, let us denote by /, g the densities of fi, u. Then, take the 

'interestingly, Knothe used this rearrangement to prove the isoperimetric inequahty, for which 
it is well suited. . .but in fact not as much as Brenier's map, which Figalli, Maggi, and PrateUi [13] 
used more recently to prove sharp isoperimetric inequahties. 
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first marginals, which we denote by /x^and v^; their respective densities are 

f^{xi) = j f{xi,X2)dx2 and g^iyi) = j g{yi,y2)dy2. 

Define as the optimal transport map between fi^ and u^. Next, consider the 
disintegration of fj, and v with respect to /J,^ and i^^, that is to say the two family 
of probabilities measures {/x^^} and {fy^} such that the densities of /x^^ and 
are 

2 , . f{xi,x2) 2 / N giyi,y2) 

For any xi let R^{xi, •) be the optimal transport map between fi^^ and ^^^1(3,^)- 
Then the rearrangement is R{xi,X2) ■= {R^ {xi) , R^ {xi, X2)). It is not difficult 
to check that it sends fj, onto v. 

What Carlier, Galichon and Santambrogio proved is that, if in Monge's problem 
the cost is, for instance, replaced with 

1 

ctix,y) = -^t'' ^\xk-ykf, 
^ k=i 

then, when the two measures are absolutely continuous, as t goes to 0, the 
corresponding optimal transport maps Tj converge in to the rearrangement 
R. When the initial measure /x is uniform and the final measure v is discrete, 
v = Y^aiSy^, they could also establish an ode governing the evolution of the 
Kantorovich potential ipt, at least when the first coordinates of the |/j are distinct. 

Thus, the following questions arise: in the continuous case, is it also possible 
to find a differential equation satisfied by V't? and if the answer is yes, is there 
uniqueness, that is to say, given the proper initial condition for t = 0, is V't 
the only solution to this equation? As we are going to see, the answer to both 
question is positive, at least, to discard boundary issues, on the torus. More 
precisely, we have the following: 

Theorem. Let At be the (1, i, . . . , t^~^) diagonal matrix, and /x, u two probabil- 
ity measures on = M^/Z-^ with smooth, strictly positive densities /, g. The 
optimal transport map for the cost 

1 ^ 

= 9 XI*'' ^dixk,ykf, 
^ fc=i 

d standing for the usual distance on T^, is then Tt{x) = x ~ A^^\/ipt{x), where 
the Kantorovich potential ipt is chosen so that / V't = 0. The map 1 1-^ ipt is smooth 
from (0, +00) to '^°°(T^), with At — D'^ipt > ^ at all times, and satifies 



div (/ [l - A^^D^i^t] ' (A^'VV't - A^MtA^^VVt)) = 0. 



(3) 
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Moreover, ifjt is the unique solution of (3) such that, if we write fort ^ 0, 

tptixi,. . . ,xn) ='tpj{xi) +t^pf {xi,X2) + ... + t^~'^ipf{xi,...,XN), (4) 
with 

Vxi, . . .,Xk-i, J iPti^i^ ■ ■ ■,Xk-i,z) dz = 0, 

then t H> (-04 , • • • , V'i^) i^ on [0, oo), and the Knothe-Rosenblatt rearrangement 
R = {R^,...,R^) is given by 

R^{xi, . . . ,Xk) = Xk - 5fcVo(^i> • • -^Xk). 

The first point is obtained by noticing that, at least when t stays away from 0, 
V't is the unique solution to a Monge-Ampere equation J^{At, ipt) = 0, where 

J^{A, u){x) := f{x) -g(x- A-^Vu{x)) det (/ - A-^D'^u{x)) , 

is defined on a proper subset of =^jv^ x ^^(T^), and then proving that we can 
apply the implicit function theorem. As it is well-known, the invertibility of the 
differential D^J' in [At, ipt) is equivalent to the existence and uniqueness of the 
solution to a strictly elliptic equation, so the argument is rather straightforward. 

For small times, because of the degeneracy of At, we need the decomposi- 
tion (4), which leads us to introduce another operator, namely: 

g{t, u\u\..., n^) := F (^t, ^ t^-^'') , 

defined on a good subset of [0, +oo) x ^2(^1) x 'rf^{T'^) x • • • x '^^(T^), in such 
away that {ipl, . . .,tp^) is the only . . . ,m^) such that ^(t, . . .,u^) = 0. 
Unfortunately, a loss of regularity for the solutions {v^, . . . , v^) of the equation 

Dug{t,^\...,^''){v\...,v'')=q 

prevents us from applying the implicit function theorem once more. We cir- 
cumvent this difficulty by using the smoothness of the Kantorovich potential ipt, 
which allows us to define ^ on a subset of [0, +oo) x ^°°(T^) x • • • x -^""(T^), 
so that to have an infinite source of regularity, and then use the Nash-Moser 
version of the implicit function theorem. 

We do not know if there is an equivalent result on M^. To be able to construct 
the Knothe rearrangement, compactness is required, but in this comes with a 
boundary. The problem is that the rearrangement is more easily contented with 
sets whose shapes are somewhat compatible with the axes, e.g. the square, but 
known regularity results for Brenier's map fail to apply in that kind of setting. 

Acknowledgements This work is part of a pho thesis supervised by Luigi Am- 
brosio (sNS, Pisa) and Filippo Santambrogio (Univ. Paris-Sud), whom the author 
would like to thank warmly for their advice and strong support. Financial sup- 
port is provided in part by a "Vinci" grant from the Franco-Italian University. 
Much of this paper is also the result of an extended stay in Pisa in Fall 2011, 
which was made possible thanks to the ens-sns exchange program. 
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2 General quadratic costs on the torus 



Given two probability measures fi, u on the torus = M^/Z^, we want to 
study the evolution with t of the optimal transport map for the cost 



1 ^ / \ 

k=l \i<k J 

where d xT^ ^ [0, +oo) is the usual distance on T^, the Aj : M — )■ [0, +oo) 
are smooth and such that Aj(t) = if and only if t = 0. For t > 0, this is a kind 
of quadratic cost on the torus. Notice that, more generally, we can define a cost 
given any positive-definite symmetric matrix A G =^jv ^ as follow: first, consider 
c : X ^ [0, +oo) defined by 

c{x,yf := M l^A{x - y - kf, 

where Az"^ is a convenient shorthand for {Az, z), and then take the induced map 
c : X [0, +oo). This is equivalent to changing the usual metric on 

with the one induced by A in the canonical set of coordinates, and then taking 
half the resulting squared distance as the cost. 

An interesting property of such a cost c is that in this case the so-called c- 
transform of a function u : — >■ M is strongly connected to the Legendre 
transform (for the scalar product induced by A) of x ^ 5^^^ ~ u{x), defined 
on (we then see n as a periodic function on M^). Let us recall that the c- 
transform of u is the map u'^ : ^ M defined by 

u^iy) = inf , {c{x, y) - u{x)} . 

This is interesting, because McCann [6] showed that, under suitable assumptions, 
the optimal transport map T can be written as T{x) = exp^(— VV'(a^)), for some 
function ip such that tp'^'^ = ip. A map u such that u'^^ = uis called c-concave. 

Lemma 1. A function u : — t- R is c-concave if and only if 



V : 



x I—)- ^Ax"^ — u{x) 



2' 

is convex and lower semi-continuous. Ifu is '^'^ and such that A — D^u > 0, then 
X I— >■ X — A~^'Vu{x) induces a diffeomorphism — t- T^. 

Proof. If u is c-concave, then v is convex and lower semi-continuous, for it can 
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be written as a Legendre transform: 



= Iax"^ - inf {c(x, y) - u^iy)} 

2 yeT^ 



sup \l^Ax^ - c{x,y) +u''{y)\ 
sup sup I -Ax'^ — -A{x — y — k)"^ + u'^{y) \ 



sup <^ (^j;, y) 



Ay^ - u\y) 



Conversely, if v is convex and lower semi-continuous, then it is equal to its double 
A-Legendre transform: 



Therefore, 

u{x) 



v{x) = sup < {Ax,y) — sup [{Az,y) — v{z)] 



-Ax' 



sup < {Ax,y) - sup [{Az,y) - v{z)] 



inf 



Ia{x - yf - I Ay' + sup [{Az,y) - v{z)] 



inf 



1 



A{x - yf 



inf 



A{y 



u{z) 



that is to say u{x) = u'^'^{x). 

If u is and such that A — D'u > 0, then by compactness A — D'u > el for 
some e > 0. Thus, v being convex with a super-linear growth, Vv : — t- 
is a diffeomorphism, and so is the map T : x x — A~^Vu{x). Notice that, 

if /c E Z^, then T{x + k) = T{x) + k, therefore T induces a diffeomorphism 

_^ jN^ □ 

In the next proposition, we start from the existence and uniqueness of the Kan- 
torovich potential for such a generalized cost (this comes from McCann [6]), and 
then apply the results of Caffarelli [14] to get its smoothness, in the exact same 
way as Cordero-Erausquin [5] did. More general results regarding the regularity 
of the potential, and thus, of the optimal transport map, on arbirary products of 
spheres have been recently obtained by Figalli, Kim, and McCann [15]. 

Proposition 2. Let fi and u be two probability measures on with smooth, 
strictly positive densities, and let c be the quadratic cost on x induced by a 
definite-positive symmetric matrix A. 
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Then there is a unique c-concave function ip : M with / ^ = such that 

T : ^ defined byT{x) := x - A'^Vipix) sends n onto v. 

The function ip is a Kantorovich potential, it is smooth, and the application ip : 
X I— >■ \Ax'^ — ip{x) is a smooth, strictly convex function on R^. 

The transport map T is optimal for the cost c. There is no other optimal transport 
plan but the one it induces. 

Of course in this proposition, instead of T{x) := x — A~^'Vip{x) we should 
have written T(x) = x — 7:{A~^'V'ip{x)), where vr : — ^ is the usual 
projection. 

Proof. Let us denote by the gradient for the metric induced by A. Then ac- 
cording to McCann [6], there is a Lipschitz function ip : M that is c- 
concave and such that T : x i— )■ exp^[—\/Aipix)] pushes /i forward to u. It 
is uniquely defined if the condition / 'tp(x) = is added, and moreover it 
is optimal for the Monge-Kantorovich problem. Notice that here on the torus, 
exp^l—VAipix)] = X — A~^Vip{x). 

For any x G M^, let ip{x) := ^Ax"^ - ip{x). Then T{x) = A'^Vip{x) sends 
fi onto u, seen as periodic measures on M^. Moreover, according to Lemma 1, 

is a convex function. Now, let V he a. open, convex subset of M^, and define 
U = {Vip)~^{V); then V99 sends onto A^ij\y, and both measures are still 
absolutely continuous with smooth, bounded, strictly positive densities. There- 
fore we are entitled to apply the results of Caffarelli [14], and thus we get that (p 
is strictly convex and smooth on U. As U is arbitrary, ip is strictly convex and 
smooth on M^. Thus, ip is also smooth, and T is a diffeomorphism. □ 

3 PDE satisfied for positive times 

Let // and u be two probability measures on with smooth, strictly positive 
densities / and g. According to Proposition 2, for any A G ^jv^> we have a 
smooth Kantorovich potential '■ — t- M. What can we say of the regularity 
of"^ : A^ ^'a? 

As X t-^ X — V^'^(x) sends /i onto the following Monge- Ampere equa- 
tion is satisfied: 

f{x) =g{x- A~^V^a{x)) det (/ - A'^DHj^ . 

For u e '^^(T^) such that A - D'^u > and / u = 0, we set 

J'iA, u) = f -g(id- A^^Vn) det (/ - A^'^D'^uj . 

Thanks to the characterization of c-concave functions from Lemma 1, and to 
Proposition 2, we have 

Lemma 3. For any u G '^^(T^) such that A — D^u > and J u = 0, we have 
J^{A, u) = if and only ifu = ^a- 
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We are now going to prove that we can apply the imphcit function theorem. 

In the following, for any function space X we denote with a o subscript the 
space formed by the elements of X having a zero mean value, e.g. is the 
space of all u £ such that J u = 0. 

Lemma 4. The operator J- is smooth. For any A e ifu € is such 

that A- D'^u > 0, ifv G 'rf^{T'^), then 



DuT{A, u)v = div (/ - T{A, u)) 



A-D' 



u 



-1 



Vv 



^ div (g (id - A~^Vuj [Co (^A-D 



detA 



T 

Vv 



We denote by M ^ the transposed matrix of M, and by Co(M) its comatrix, 
that is to say the matrix formed by the cofactors. 

Proof. The smoothness of T is obvious. By substitution, for any ^ G 

A-^\7u{x)) [fix) - HA, u){x)] dx = I i{y)g{y) dy. 

Therefore, if we conveniently set T^u(x) := x — A^^Vn(x) and differentiate the 
previous equation with respect to u along the direction v, we get 

- j (v^Tau), A-^Vv) if - T{A, u)) - I ^{Tau)D^T{A, u)v = 0. 

Since o Ty^u] = [DTy^ul'^V^iTj^u), we have 

(vC{Tau),A-^Vv) = (v[^oTau],[DTau]-^A-^Vv) 

= (V[C o Tau], [I - A-^D\]-^A-^Vv) , 

and this yields 

j ^{TAu)D^F{A,i^A)v 

= j ^Tau) div ((/ - J'iA, u)) [I - A-'D\]-^A-^Vv) , 

and thus, since ^ o Tau is arbitrary, we get the first equality. Then, we can easily 
obtain the second expression using the formula = [Co M]^/ det(M). □ 

Lemma 5. Lete > and A € ■ Ifu G ^^^(T^) is such that 

A-D^u > e{detA)^I, 
then for any q G [Hl{T^)]*, there is a unique v G Hl{T^) such that 

DuT{A, u)v = q. (6) 
Moreover, \\v\\ifi < Cc\\q\\(^Hi-^*. 
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Proof. As A- D'^u > e{det A)^/^^'^h, the lowest eigenvalue of Co {A - D'^u) 
is bounded by e^^^ det A. Since g > 5 ior some (5 > 0, for any ^ G '^°°(T^), 



det A j |V^|2 < y"([Co(^-Z)2n)]TvC,V0 

]j9{^d- A-'Vu) ([Co - I)2y)]^Ve, VO, 



< 
- 5 



and thus 



/ |VeP < / ^D^HA,u)(. (7) 

Therefore, thanks to the existence of a Poincare inequality on H^(T^), the map 
iC^v) ^ I '<lDu^{A,u)^ induces a coercive, continuous bilinear form on H^. 
We are thus entitled to apply the Lax-Milgram theorem, which yields the ex- 
istence and the uniqueness, for every q G (Hi)*, of a t; G satisfying (6). 
Moreover, (7) immediately gives us < ^jv-i lkll(_H'i)*- D 

The regularity of the solutions to an elliptic equation is well known. However, 
as in the following we will need some very precise estimates to apply the Nash- 
Moser theorem, let us give a proof of the following result: 

Lemma 6. Under the same assumptions, and with the same notations, for any 

n>l, ifu G "^"^^ andq G H^~^ satisfy \\u\\^g3 + \\q\\(^Hl)* < thenv G 
and 

\\v\\jjn+l < C^^M.n + ll'"ll<!f"+2} • (8) 

Proof. We proceed by induction. Let n > 1, u G "^"^^ and q G such that 
A - D'^u > e(det and ||u||c^3 + < M. We assume that we 

already know that the corresponding solution v is in H^, and that 

IHIh" < Ce,M,n-l{\\q\\H^-2 + \\u\yn+i} . (9) 

Notice that we do have such an inequality for n = 1, according to the previous 
lemma, but with instead of Let us now show that it implies 

V G i^^+i and 



First, we set B^u := (/ — J- {A, u))[A — D^u]^^, so that Equation (6) becomes 

DuJ'iA, u)v = diY{BAuVv). (10) 
Then, for /i G and ^ G we also define 

Th^{x) := ^{x + h) and Oht.{x) := . 

n 

Notice then that 6h{vO = V^h^ + (4??)^/,^ and II^^IIl^ < U\\h^. 
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Let g be a 2-index, with \v\ := vi + V2 = n — 1, and let /i £ be small 
enough. We can apply the operator 5h to Equation (10), and we then obtain 

(i\\{BAuV5hv) = 5hq - div [{5hBAu)VThv] 

Then, by applying d,^, we get 

div {BAuVShd^v) = Shd^q - | j div [{5hdv-aBAu) Vrhdav] . 

- E \]<^^^[{du-aBAu)V5hd^v\. (11) 
Now, Lemma 5 tells us that this implies 

\\^hduV\\Hi <Ce\\5hdvq\\(^Hl)* 

+ Ce ^ I j ||div [{5hdy-a.BAu) ^ThdaVlW^^^ly 

+ Ce I j ||div [{dy-aBAu) y^hdaVlW^^Hiy . 

0<a<u V"/ 

Since < ||5jy(?||L2, this bound is uniform in h, and so it is enough 

to ensure v £ H^^^ and 

||w||_H-"+i < C* < + ^ (1 + lkll'^"-fe+2)ll^ll_fffe+i [• • (12) 

[ 0<A;<n-l J 

Notice that, when n > 1, the following Landau-Kolmogorov inequalities hold 

||tf||^n-fc+2 < C'fc_„||u||<:^3 " ||u||f^„+2) 
k ^ _ k 

\\v\\^k + l < Cfc^n||f llj^i llflljlfn" 

They are quite classical and can be easily proved by induction from 

mv^^ < and n^n^i < ^uhmH^, 

for ^ smooth enough satisfying / C = 0. Since a^^^h'' < (1 — t)a + tb, we get 

,, ,, ,, ,, k / k \ 

11 — 1 \ n — 1 J 

and therefore 
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This last inequality still holds when n = 1, thanks to (12). In any case, as 
||f < Celkll/f-i and ||n||<:^3 + < M, using our assumption (9), 

||t;||j^n+i < Cs,M,n{\\q\\H"-'^ + • 
This is exactly what we wanted. □ 

Lemma 7. Under the same assumptions, for any q G '^'"(T^), there is a unique 

V e <^<r+^'"(T^) such that 

DuF{A,u)v = q. 

Proof. If g G 'rfo'°', then q £ H^, and thus according to the previous lemmas, 
there is t; G such that DuJ^{A, u)v = q in {Hi)*. But since J q = 0, given 

the particular form of DuJ^{A, u)v given by Lemma 4, such an equality in fact 
holds in H~^. Thus, locally, in a weak sense, 

DuJ^{A,u)v = q. 

Then, we can locally use the theory of regularity for the solutions to a strictly 
elliptic equation in to get existence and uniqueness of u G (cf for 

instance Gilbarg & Trudinger [16], Chapter 6). □ 

Theorem 8. For any A G ^jy^' Kantorovich potential between 

the probability measure fi and v, which are still assumed to have smooth, strictly 
positive densities. Then, for anyn > and a G (0, 1), the following map 

Proof We denote by n be the set of all {A, u) G x <r"+^'°(T^) such that 

A - D^u > 0. Then n is open, the operator J" : S7 ^ defined by 

T{A, u) = f -g{id- A-^'Vuj det (/ - A^'^D'^u 

is smooth and, according Lemma 7, DuJ'{A, t/^A) '■ '^""'"^'"(TT^) — t- ^^'"(T^) is 
a bijection. From the Banach-Schauder theorem, we deduce it is an isomorphism. 
Since J- {A, ^'^) = 0, according to the implicit function theorem, there is a 
map $ defined in a neighborhood of A such that B — D^^b > and, for any 
u G '^<r^^'", J-{B,u) = if and only if u = (^b- According to Lemma 3, it 
impHes ^b = ^b- Thus, globally, ^' = ^> is a map "^"'"(T^). □ 

We are now going to apply this result to the cost c defined by (5), that is to say 
the cost induced by the matrix 

/ 1 \ 

Ai(t) 

At := Ai(t)A2(t) 

V J 
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where Ai, . . . , A^v-i : M — >■ [0, +00) are assumed to be such that Afc(t) = if 
and only if t = 0. 

Theorem 9. IfXi, . . . , Xn-i are smooth, the mapip : t ^j^^ is^^, and satisfies: 



div I / U - D^i,t] (vV't - AtAj^V^t) ] = 



(13) 



Moreover, ifu : (0, +00) c^n+2,a^jN^ j^^i and satisfies, for all t G (0, +00), 



At-D'^u > and div < f At - D'^ut (Viit - AtA^^Vut 



-1 



= 0, (14) 

and ifutQ = ipto for some to > 0, then ut = ipt far all t > 0. 

Proof. If ipt '■= ^At^ foi' all t > 0, we have J- {At, ipt) = 0. If we differentiate with 
respect to t, we get 

DuF{At, iJt)i^t + DAF{At, iJt)At = 0. 
We have seen in Lemma 4 that 



DuJ'{At,i;t)i^t = div ( / At- D^^Jt 
On the other hand. 



DAF(At, iPt)At = - div (j [At - D^^Pt] ' AtA-'ViPt 
We thus get (13). 

Ifu : (0, +00) ^n+2,a^jN-^ ^1 ^^^j satisfies (14), with uto = V'io for some 
to > 0, then J^{At,ut) must be constant and equal to J^^ATQjUtg) = 0. Thus, 
according to Lemma 3, ut = "^At ■ 



4 Initial condition in dimension 2 

Due to the very technical nature of the proofs, the following sections will only 
deal with the dimension 2, to ease the exposition. In the last section, we shall 
explain what changes in higher dimension. 

Let A : M — 7- [0, +oo) be a smooth function such that A^ = if and only if 
t = 0. From now on, we will only consider the cost induced by 



At 



1 
Xt 



For t 7^ 0, let ipt be the associated Kantorovich potential between the probability 
measures fi, v, assuming they have the same properties as before (that is, strictly 
positive and smooth densities), and let Tt be the corresponding optimal transport 
map. Then, according to Theorem 9, t ^ ipt and t ^ Tt are on R \ {0}. 
Moreover, Carlier, Galichon, and Santambrogio [10] proved: 
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Theorem 10 (C.-G.-S.). Ast ^ 0, the mapTt converges to the Knothe-Rosenblatt 
rearrangement R in L^(T^, /x; T^). 

Let us denote by Uq{-) and Uq{xi, •) the Kantorovich potentials for respectively 
-R^(-) and •). Indeed, recall that sends the first part fi^ of the disinte- 

gration of /i onto the first part of the disintegration of z/, and that R^(xi, •) 
sends the second part onto i^^ij-^.^), in an optimal way for the squared dis- 
tance on the 1 -dimensional torus T^; hence these transport maps come from 
some potentials. We have: 



R{x) 



xi - diuKxi) 

X2 - d2ul{xi,X2) 



The Carlier-Galichon-Santambrogio theorem suggests some connexion exists 
between -0* and {uq,Uq). Since = id — A'^^Vijjt, let us follow our instinct 
and set 

'lpt{xi,X2) = Iplixi) + XtTpt{xi,X2), 

and, to ensure uniqueness, require 

J V't(xi)dxi = and J ^f{xi,X2)dx2 = 0. 
Notice that tpj and ipf are then uniquely determined, and are smooth, since 

= /«....)d.. and =i 
Let us denote by E the set of all (t, u\ u^) G M x '^'^(T^) x <^°°(T2) such that 
/n^(x,)d.,=0 and / d.. = 0, 

and by the open subset of E formed by the {t,u^,u'^) such that: 

• either t / 0, and then At - D'^{u^ + XtV?) > 0; 

• or t = 0, and then 1 — di^iu^ > and 1 — d2^2U^ > 0. 

Then, thanks to Lemma 1, we can define an operator Q : ^} ^ '^°°(T^) by 
setting, when i / 0, 

g{t, u\u'^) := J^{At,u^ + Xtu^), (15) 
where is the operator introduced in Section 3: 

J^{A, u) = f -g{id- A-^Vuj det (/ - A'^D'^uj . 
Since, according to Lemma 3, J-'{At, u) = if and only if n = ipt, we have: 
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Lemma 11. For any (t, u^, u^) G G{t, u^, v?) = if and only ifu^ = -0^ and 
Now, we are going to extend Q for t = 0. Notice indeed that 

A V(n +Xtu)=[^ ^^^2 j 

and A D {u +X,u)=^ ' j , 

therefore we can smoothly extend Q. If we conveniently define an operator d by 
setting du := {diu^, d2V?), then = id — OiIjq, and 

g(0, u^,u^) = f -g (id - du) det (/ - Ddu) . (16) 

Alas, we cannot do the same as in the previous section and apply the implicit 
function theorem, for if we solve D„^(0, t/'q, V'o)(^^) = 1^ then a priori the 
solution is not smooth enough. Indeed, as we will see later, if g G H^, then 
^1 g \)^^ -yye can only get u ^ G i/". However, we can circumvent this dif- 

ficulty by considering functions, so as to have an infinite source of smooth- 
ness, and use the Nash-Moser implicit function theorem instead of the "classical" 
implicit function theorem. 

Before stating our next result, let us recall some definitions from the Nash- 
Moser theory. For more details, see for instance Hamilton [17]. 

Let X and Y be two Frechet spaces, endowed each one with a family of in- 
creasingly stronger semi-norms, denoted by {|| • ||i^}n>o and {|| • ||^}n>o- For 
instance, you can think of ^°^(T^), with the norms || • ||„ = || • ||<:^n or equiva- 
lently || • ||„ = || • \\h^. A map 99 : [/ — )■ y is said to be "tame" if it is defined on 
an open set U <Z X,is continuous, and in a neighborhood V of each point, one 
can find r > 0, 6 > and a sequence {Cn)n>b of positive constants such that the 
following "tame estimate" is satisfied: 

Vx G y,Vn > h, Mx)\\l <Cn{l + . 

Notice that r, h, Cn can depend on V, but V cannot change with n. The map (p is 
"smooth tame" if it is smooth and if all its Gateaux derivative D^(p : [/ x X — t- y 
are tame. From the definition (15) of Q, we easily get: 

Lemma 12. The operator Q : 17 ^°°(T^) is smooth tame. 

The Nash-Moser theorem holds for some Frechet spaces, the so-called "tame 
spaces" defined as follows. If is a Banach space, the space of exponentially 
decreasing sequences in E is defined as: 

E(i?) := j(u„) G i?^ ; Vn G N, := £ e^''\\uk\f < ooj 

I k=0 J 
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A Frechet space X is said to be "tame" is there is a Banach space E and two tame 
linear maps $ : X !;(£') and : T^{E) X such that o ^ = idx- For 
instance, ^°°(T^) is a tame space. If X and Y are tame, then so is their cartesian 
product X xY. 

Theorem 13 (Nash-Moser). Let X and Y be two tame spaces. Let U d X he an 
open subset and ^ : U Y be a smooth tame map. We assume that, for any u & U 
and any q G Y, there is a unique v G X such that D^{u)v = q. If the inverse 
operator S : U x Y X is a smooth tame map, then $ is locally invertible, and 
the local inverse maps are smooth tame. 

Corollary 14 (imphcit function). Let X, Y, Z be three tame spaces, and U C X 
and V C Y be open subsets. We assume ^ : U x V ^ Z is a smooth tame map 
such that ^{uo,vo) = for some {uo,vo) E U xV. If, for any u E U, v E V 
andq G Z, there is a unique w & Y such that Dy^(u,v)w = q, and if the inverse 
operators : UxVxZ^Yisa smooth tame map, then there is a smooth tame 
map defined in a neighborhood ofuQ and taking values in a neighborhood ofvQ 
such that $(it, v) = if and only ifv = ipiu). 

Notice that we need only to use this last statement in an open neighborhood 
of (0 , Uq, Uq) G ri, with Uq and Uq the Kantorovich potential associated with the 
Knothe-Rosenblatt rearrangement. Let us define this neighborhood Qq in the 
following way: First, take e > such that 1 — di^iUQ > e and 1 — d2,2Uo > 
Then we take for fio the set of all (t, , u|) G such that: 

ift = 0, l-di^iv}>e and 1 - ^2,2"^ > e, (17) 

and 

if i 7^0, 1 - di^iv} - Xtdi^iu^ > £ and A* - ^^(u^ + Aju^) > eA*. (18) 

Recall that we denote with a o subscript the sets of maps with zero mean value: 
'io^ is thus the set formed by the smooth functions u such that / u = 0. When X 
is a 2-variable function space, we also denote by a "*, o" subscript, as in ^h.^CT^) 
the set formed by the ^ G X such that / ^(-j X2) dx2 = 0. 

Theorem 15. For all {t,u^,u^) G Oq, for any q G ^^°°(T^), there is a unique 

{v^,v^) G '^^{T^) X "^,^(¥2) such that 

Dugit,u\u^){v\v^)=q, (19) 
Moreover, the inverse operator 

f nox ^-(T2) ^ ^-(Ti) X <r-(T2) 

is smooth tame. 
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Proof. We will show the existence of {v^,v'^) in Section 5. We also report the 
proof of the existence of a tame estimate for the inverse operator S to Section 6. 

Let us conclude from that point. Now all that remains to show is that S is 
continuous, and that the derivative D^S are tame. 

First, if (tfc, M^, g^) € l^o converges towards (t, u^^u^, q) € i7o> foi' each k 
let {v\,vf.) be the corresponding inverse. Thanks to the tame estimate from Sec- 
tion 6, v\ and vf, are bounded in all the spaces Hence, compact embeddings 
provide convergence, up to an extraction, to some v^,v^ as strongly as we want, 
which, since DQ is continuous, must be the solution of DQ{t, u^, u^)(w^, v^) = q. 

Then, all the derivative D^'S are also tame, since they give the solution to the 
same kind of equation as (19). Indeed, by differentiating (19), we get 

D^gDS = Dq-D{D^g), 

and then we can apply the results of Section 6 once more. □ 

If we now set tpQ = Uq and ipQ = Uq, with Uq and Uq the Kantorovich potentials 
for the Knothe-Rosenblatt rearrangement, we can state the following: 

, '^i 2\ smooth. 

Proof. On some interval (— r, r), this is a direct consequence of Corollary 14, 
Theorem 15, and Lemma 11. For larger t, it follows from Theorem 8. □ 

Theorem 17. The curve formed by the Kantorovich potentials {ipt) is the only curve 
in "^^{T^) defined on R such that, fort / 0, 

At - D'^^l^t > and div (f U - D^^JtV^ fv^t - AtA^^Vt^t)] = 0, (20) 



and that can be decomposed into two smooth curves (ipt) '^^'^ ("0?) ^uch that 

ll)t{xi,X2) = 1pl{xi) + AtV-t (X1,X2), 

with ipQ andipQ being the Kantorovich potentials for the Knothe rearrangement. 

Proof. Let Ut = + XfU^ be such a curve, and let us check that ut = ipt- Since 
Uq and Uq are the potentials for the Knothe rearrangement, (0, Uq, Uq) G ^q, so 
{t, ul,uf) is in J7o least for t small. For t ^ 0, (20) is equivalent to 

DuF{t, ut)ut + DtF{t, Ut) = 0, 

and therefore 

DuG{t,u\,ul){u],ul) + DtG{t,u\,ul) = 0. 

By assumption, ^(0,Uq,Uq) = 0. Integrating in time, we get Q{t,u^ ju"^) = 0. 
Therefore, according to Lemma 11, = tpt ^^'^ ''^t — ''Pt' i-^- '^t = i^t- 
For larger t, we apply Theorem 9. □ 
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5 Proof of the invertibility 

We recall that 



F{A,u) = f -gild- A-'Vu] det / - A-'D^u] , 



and 



with 



At :-- 



1 

Xt 



(21) 



We want to prove the invertibility of DuG{t, n^, n-^). The first lemma will con- 
sider the case t ^ 0, the second one the case t = 0. 

Lemma 18. For any (t, u^, u^) G i7o with t ^ 0, for any q S ^^(T^), there is a 
unique {v\v'^) G <^^°°(T1) x '^^'^(T^) such that 

Dug{t,u\u^){v\v^) =q. (22) 

Proof. Let {t,u^,u'^) G l^o with t / 0, and let q G 'tf^{T^). Then, if we set 
:= + AtU^, Lemma 6 tells us that there is a unique G '^<f°(T^) such that 

-1 



div /-g(t,n\n^) I - At'D^ut A^^Vvt 



(23) 



Let us define 



Then, by construction, {v^, v'^) is the unique pair solving (22). 



□ 



Lemma 19. For any {0,u^,u'^) G ^o> for any q G '^<^(T^), t/iere is a unique 

{v'^,v^) G "tf^iT^) X '^^'^(T^) swc/i f/iflf 

Dug{0,u\u^){v\v^) =q. 

Proof. We want to solve 

D^g{0,u\u^){v\v^) =q. 
By substitution, for any ^ G Equation (16) yields 

ax-du{x)) [f{x)-g{0,u\u^){x)] dx = ji{y)g{y)dy. 

Therefore, if we differentiate the previous equation with respect to u along the 
direction v, and recall our notation du = (5in^, d2V?) and dv = (9if ^, d2v'^), 
we get 

(Ve(id - du),dv) (/ - g{0, u\u'^)) 

^(id - du)Dug{0, u^,u^){v^,v'^) = 0. 
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Since o (id — du)] = [/ — Ddu\'^V^{id — du), we have 

{Vi{id-du),dv) = (v[io {id - du)],[I - Ddu]-^dv'^ 
and this yields 

DuQ{Q,u^,u^){v^,v'^) = div((/-g(0,u\u2)) [I ~~ Ddu\-^ dv) . 
Notice then that 

(/-e(o,«'.„'))[/-o8„l-' = s(id-8u)( i-^<^^f i-a°,,„' )■ 

thus, 

= di g {x - du{x)) - d2,2u'^{x)j div^{xi) + ^2 [. . .] . 
Therefore, if DuQ{0, u^,u'^){v^, v"^) = q, integrating with respect to X2 yields 



di g{x — du{x)) — d2,2u^{x)j div^{xi) 
that is to say 



dxo 



q{x) dx2, 



di 



' {x — du{x)) (^1 — d2,2u'^{x)J dx2^ div^{xi) = Jq{x)dx2- (24) 

But there is a smooth map Q : — )■ M such that diQ{xi) = J q{xi, X2) dx2, 
since J q{x) dx = 0, and it is unique if we require Q{0) = 0. Thus, taking a 
primitive of (24), there is a c G M such that: 

' (x - du{x)) (1 - d2,2U^{x)j dx2 div^{xi) = Q{xi) + c. 



Since G{xi) > 0, we get 



Q + c 
G ' 



and this yields the unique possible value for c since the integral w.r.t. xi of the 
right hand side must be zero. Combined with the condition / dxi = 0, we 
thus have completely characterized v^. 
Now, let us do the same for v'^. We have to solve the equation 



d2 



g{id-du) [l-di^iu^j d2V^ 
- q — di 5 (id — du) ( 1 — 82^2^^ ) div^ 



do 



g (id — du) di^2U^div^ 



and this is exactly the same kind of equation as (24). If we fix xi G T^, the same 
reasoning can be applied here, and thus we get v"^. □ 

This ends the proof of the invertibility. All that is left to show is that we have 
some tame estimates. 
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6 Proof of the tame estimates 



Our aim here is to show that, locally on {t,u^,u'^) G and q G '^^(T^), for 
any n E N, there is a constant C„ > such that, if 

Dug{t,u\u^){v\v^) =q (25) 

for some {v^,v^) G ^^°°(Ti) x ^.'^(T^), then 

||u^||//n + 2 + IIV^II//" < Cn(l + \t\ + ||u"'"||/fn + 3 + Hu^H/fn + S + HgH/f") • 

In fact, we will prove something slightly stronger: 

||f"^||j:/n+2 + ||92V^||_ff" < Cn (||'U"'"||'^"+3 + ||n^||<g'n+3 + HgH//") • (26) 

Indeed, since J v'^{xi, X2) dx2 = 0, we have a Poincare inequality, which im- 
plies < Cn||02t'^||H"- Notice also that (26) would by itself prove there is 
uniqueness. 

We start with the case t 7^ 0. As the bound for ||f ^||j:^n+2 simply follows from 
Lemma 6 and an integration with respect to X2, we just have to find a bound for 
||52f ^ll/f". Let us begin with ||92f^||i2. 

Lemma 20. Let M,e > 0. There are C = C{M, e) such that, if{t, u^,u^) G 
andq G '^^{T^) satisfy 

\\q\\L2 + Wu^Wc^s + llii^ll^fs < M, (27) 

if{v\v'^) G ^^°°(Ti) X <^^°°(T2) is a solution of (25), then 

\\d2V^\\L2 < C. (28) 

Proof. We set Ut := + Xtu^ and also vt := v"^ + Xtv"^. Then, by assumption, 

DuF{At,ut)vt = q. 

The property (18) in the definition of 0,q ensures we an apply Lemma 6 and get 



\vt\\m < C'e,A',l{lkllL2 + ||Mt||<r3} < C. 



(29) 



We now set 



Bt:= if -g{t,u\u^)) I-A^'D'ut 



1-1 n2„ 



-1 



A 



-1 



f-g{t,u\u') 

det{At - D^ut) 
g(id - At^Vut) 
det At 



Co {At - D^ut) 



Co {At - D\t) 
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so that, according to (21) and Lemma 4, Equation (25) becomes 

div{BtVvt) = q. 

Notice that det At = Xt and 



Co {At - D\t) = { ^* ~t'^''f .^'^T^ I ' 



therefore we can write 

Bt = Ut + Vt/\t (30) 



with 



Thus, 



f/.:=,(id-A-V..)(^-^^^f Y j' 

y.:=,(id-A-V..)(0 (32) 



As ^ = 0, we have 14Vw^ = 0. Since vt = v"^ + Xtv'^, we get 

diy{UtVvt) + diy{VtVv^) = q, 

that is to say 



do 



g{id - A^^Vut){l - di,iut)d2v^] =q- div{UtVvt). (33) 



Since g > S for some 5, and as (18) in the definition of means 1 — di^iut > £, 
allowing the constant C to change from line to line we get 



,.2|2 



\\d2vYL^ < ^ / 5(id - A-'Vut){l - di,iut)\d2v' 
<C J[q- div{UtVvt)] 
<C{\\qh2 + \\UtVvt\\H^)\\v''h2 
However, since / v'^{xi,X2) dx2 = 0, we have ||w^||l2 < C||92f^||i2. Therefore, 

\\d2V^L2\\v^L2 < C\\d2vYL2 < C{\\q\\L'2 + \\UtVvt\\H^) \\v^l^. 

Thus, since ||f/t||f^i < C(l + ||^i"'"||<^^3 + H'li^ll'g's) < C as we can see from (31), 

\\d2v'\\L2<C{\\q\\L2 + \\vt\\H2}. 

Then, using (29), we get the result. □ 
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We now proceed by induction to get an estimate for any order n G N. 



Lemma 21. Under the same assumptions than in the previous lemma, for any n G 
N, there is a constant Cn = Cn{M, e) such that 



\\d2v'^\\wr < Cn (ikll//" + ll'""'"ll<r"+3 + ll^^lh^^+a) 



(34) 



Proof. Let us assume (34) has been proved for some n G N, and let us show it 
holds even for n + 1. Let G be such that |z^| := ui -\- 1^2 = n -\- 1. Recall (33), 
that is to say 



do 



g{id - A-^Vut){l - di,iut)d2v^] =q- dw{UtVvt). 



We already know from Lemma 7 that vt = + Xtv"^ is smooth, therefore, if we 
apply du, we get 



do 



82 



E 

0<a<v 



+ dyq - du div(C/tVft) 
On the other hand, since g > 5 and 1 — di^\Ut > e, we have 

\\d2duvYv^ < ^ / 9(id - A~^Vut){l - di,iut)\d2dy\^ 

'g{id - A-'Vut){l - di,iut)d2d,v^^ 



< 



6e 



do 



duv\ 



Thus, 

\\d2dyV^\\\2 

< E I 

0<o<i^ 

and therefore 



d2dyv'^ 



[dyq - du dw{UtVvt)] dyV^ 



1 

Te 



\d2duvYL^ 

< c\\du-a{g{id-A-^Vut){l-di,iut)]d2da 



0<a<u 



L2 



d2dvV' 



C\\duq-dudiw{UtVvt)\\L2 



dy 



L2 



L2 
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As ||9,yV^||j;^2 < c||529;^w^||l2, we get 



\d2duV^\\L2<C l5(id- A-iVnt)(l-9i,int) 

0<fc<n 



c^n + l — fc 



+ C{\\q\\H^+i + \\UtVvt\\H^+2] . (35) 

On the one hand, we can use the same Landau-Kolmogorov inequalities as in the 
proof of Lemma 6, and use again the fact that a^^*6* < (1 — t)a + tb, to get, for 
< < n, the following bound: 



<cn{\g{id- A~^Vut){l-di,i 



L2 



g{id- A-'Vut){l-di^iUt) d2v' 



Recall we have assumed (34) holds true for n, therefore, using (27), we get 



g{id - A-^Vut){l - di,iut] 



c^n + l — k 



d2V^ 



JJk 



(36) 



On the other hand. 



< c {\\Utyn+2\\vvt\\Hi + \\Ut\yi\\vvt\\H^+2} , 

which, since ||^i^||c^3 + ||^t^||'g'3 < M, implies 

\\UtVVt\\H" + 2 < C|(^l + ||u"'"||<^n+4 + ||ti^||%fn+4j 1 1 ^'t 1 1 //2 + ||Ut||^n + 2| . 

Then, using Lemma 6 we get 

\\UtVvt\\H"+2 < c„ (II^IIh" + + lk^ll<r"+4) • (37) 

Bringing together (35), (36), and (37), we get the estimate we sought. □ 
Lemma 22. The result of Lemma 21 still stands fort = 0, with the same constants. 
Proof Let (0, u^,u'^) e and q G "tf^iT^) such that 

lkllL2 + h^Wv'-i + \\u^\W3 < M, (38) 

Then, since (s, u^, u^) £ ^Iq for s small enough, we can proceed by approxima- 
tion. Indeed, if {v].,Vg) is the solution to 

DuQ{s,u^,u^){vl,v'^) = q, 
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where u^^u^^q have been all fixed, then all the norms of vl , v'^ are bounded 
according to Lemma 21. Up to an extraction, there is convergence, which by 
compact embedding is as strong as we want. But the convergence can only be 
towards the solution of 

D^g{0,u\u^){v\v^) =q, 
hence estimate (34) is still valid for the limit. □ 

This proves the existence of tame estimates, at least in dimension 2. Let us now 
see what changes in higher dimension. 

7 Higher dimension 

The difficulty in extending those results in higher dimension only comes from 
the technical nature of Sections 5 and 6. We need a decomposition, not only of 
the potential, but also of the field matrix B, extending (30). The existence of such 
a decomposition is the only new difficulty. 

Setting and notations 

We consider Ai, . . . , Aat-i : M — )■ [0, +oo), assumed to be smooth and such that 
Afc = if and only if t = 0. We then define At by 



/ 1 



A:-- 



\ 



Ai(t) 



Ai(t)A2(t) 



V um J 

The decomposition of the Kantorovich potential ipt becomes 

i/jtixi, . . .,xn) = ipK^i) + AiV'(^(a;i,a;2) + . . . + Y[ >^i] ipf^{xi, ■ ■ ■,xn) 

\i<N ) 

where -(/^^ depends only on the k first variables xi, . . . , x^, and is such that 



Vxi, . . . ,Xfc_i, 
For convenience, we set 



i)\{xx, ■ ■ .,Xk-i,yk) dyk = 0. 



\k<i<N 



N 
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so that we have 
and 

Vxi,...,a;fc_i, j "' j '^ti^i-^ ■ ■ ■ ^^k-uVk, ■ ■ ■ ,yN) dyk . . . dyN = 0. 
For instance, if d = 3, we have 



= + AiV-l + AiAaV-l and 



t/jl = iPl + AiV'/ + A1A2V'? 



V'i^ = V/ + A2Vf 



Let us denote by E the set of all (t, u\ . . . , -u^) G M x n '^""(T'^) such that 

VA; G {l,...,iV}, Ju'' dxk = 0. 

Then, if {t,u^, ... , u^) G E, we set 

u^:=u^, 2^= := tx*^ + Afca'=+\ u := u\ 
and this is consistent with the previous notation. Notice that 



Vu 



and thus. 



and A ^Vu = du 



( diov} 



A-^D^u = Ddu 







( diu^ \ 

d2U^ 
V dNU^ J 

\ 




(39) 



y di^NU^ d2,NU^ ■ ■ ■ dN,NU^ J 

We define O as the open subset of E formed by the {t, u) such that: 

• either t 7^ 0, and then At - D'^u > 0; 

• or f = 0, and then 1 — dk,ku'^ > for all k. 

As previously, we need only to work on a neighborhood f^o of the Kantorovich 
potential (0, Uq, Uq), which we will define precisely later. 
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Invertibility 

We want to solve, for (0, u) € 0,q, the equation DuG{0, u)v = q. Since for t > 0, 

-1 



I -A-^D'^u 



DuQ{t,u)v = diw ^{f-g{t,u)) 

which, when replacing A~^D'^u and A'^Vv with Ddu and dv, becomes 

DuGit, u)v = div ((/ - g{t, u)) [I - Ddu]-^ dd] 



and since, when t = 0, we have 
solve is 



and du = du, what we would like to 



q = DuGiO, u)v = div ((/ - g{0, u)) [I - Ddu]~^ dv) . 



The trick is to integrate with respect to Xk+i, . . . ,xi\i to get an equation on 



, v". If 



k-1 



have already been found, [/ — Ddu]^ being lower 
triangular thanks to (39), the resulting equation on is of the same kind as the 
one we have dealt with in Section 5. The same reasoning can thus be applied. 



Tame estimate 

As in the 2-dimensional case, we need only to find a tame estimate when t ^ 
for the solution {v^, . . . , v^) of 



q = div{B\7v) with B 



g{id - A-^Vu) 
det^ 



Co {A - 



First, notice that by integrating with respect to xj^, we are reduced to the — 1 
dimensional case. Therefore, we can proceed by induction on N. 

So let us assume we already have a tame estimate for v^, . . . , v^^^. To get an 
estimate for = v'^ , we will find one for each v'', this time by induction on 
k. Since = v satisfies a nice strictly elliptic equation, and thus comes with a 
tame estimate, we need only to show how to get one for -u^ if we have one for 



^1 ^^—1 
V . . ,v . 



The key lies in the following decomposition of the matrix B: for any k, 

1 1 



B = B^ + ^B^ + -^B^ + ... + - 

Ai A1A2 Ai • • • Afc_2 



At--A 



B'', 



k~i 



where the coefficients (&^^) of B'^ are zero except when min(a,/3) = i, and 
where the coefficients (6^ g) of i?'^ are zero except for min(a, P) > k : 



( 



B' 



\ 



( 



\ 



B^ = 




\k ■ 


■ \n 
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the point being that all the coefficients b^^t p^^a /3 '■'^^ bounded in by the 
norms of the in uniformly in t, at least for small t. Indeed, if such a 

decomposition exists, since 

v = v^ + Aiw^ + . . . + Ai • • • Xi-2v'~'^ + Ai • • • \i^iv\ 

with diV^ = if i > j, which implies diV = Ai • • • Xi-idiV^, we have 



diY{BVv) 
and thus 



E x \ div(j?-V^) 



+ ^ div(S^Vu), 

Ai • • • Afe_i 



div{BVv) 



i<k 



(40) 



On the one hand, the matrix B^ is symmetric and non-negative, and we define 
r^o such that 



On the other hand, since 
Vxi, . . . ,Xk~i, 
we have a Poincare inequality: 



. . . , Xjq) dXk ... dXAT = 0, 



i>k 



2 

L2 



Therefore, 



< ^ / {&Vv^,Vv'') < ^ ||div(5^Vu^; 



2 

i2 - e 



L2 



L2 ' 



and this shows how we can deduce a estimate for from (40) and a series of 
estimates for v^,i < k. Estimates for the norms H^, n > 0, easily follow, by the 
same reasoning as in Section 6. 
Thus, all we need is the decomposition 



1 



1 



with 



B = B' + —B' + + . . . + 

Ai A1A2 

/ 



1 



B' 



Ai • • • Aat-i 



B^, 



-'i,N 
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Remember that 



B := 



g{id - A-^Vu) 
detA 



Co {A - D'^u) 



and det A = \i (A1A2) • • • (Ai • • • therefore all we have to do is to show 

how in Co {A — D'^u) we can gather the A^ so as to get the decomposition we 
seek. Since diju = Ai • • • Xras;>,{i,j)-idi,jU^'"'^'"^\ 



Co {A - D'^u)] ..= J2 n - ^^") 



k,a{k) 



(tG6„ l<k<N 
cr{i)=j kj^i 



Yl ^1" ' -^max(fe,cr(fe))-l 

creGn l<k<N 
a{i)=j ky^i 



Thus, if i < j, we set w^,^ = Aq • • • Amax(a,^)-i (^a,/3 - da,pu'^^''^°''^^y and 
then we get 



CoiA-D^u)], ,= J2 n >^i---h-i^k,a 



(fc) 



= E 

that is to say, 

Co {A - D^u) 



cr66„ l<A:<iV 
a{i)=j k^i 



Ai • • • Xi-i 



Ai---Afc_i 

i<fe<Ar 



n ^k,a{k) 
l<k<N 
L fc^i 



del A 
Ai • • • Xi-i 



J2 ^(^) n ^k,a{k)- 



o-66„ 



1<A;<A' 



(T{i)=j ky^i 

Since we have assume i < j, this is exactly what we wanted. 



References 

[1] G. MoNGE, "Memoire sur la theorie des deblais et des remblais", Histoire de 
I'Academie royale des sciences de Paris, pp. 666-704, 1784. 

[2] L. Kantorovich, "On the translocation of masses". Journal of Mathematical 
Sciences, vol. 133, pp. 1381-1382, 2006. The original paper was published in 
Dokl. Akad. Nauk SSSR, vol. 37, no. 7-8, pp. 227-229, 1942. 

[3] Y. Brenier, "Decomposition polaire et rearrangement monotone des 
champs de vecteurs", C. R. Acad. Sci. Paris Ser. I Math., vol. 305, no. 19, 
pp. 805-808, 1987. 



28 



[4] , "Polar factorization and monotone rearrangement of vector-valued 

functions", Comm. PureAppl. Math., vol. 44, no. 4, pp. 375-417, 1991. 

[5] D. Cordero-Erausquin, "Sur le transport de mesures periodiques", C. R. 
Acad. Sci. Paris Ser. I Math., vol. 329, no. 3, pp. 199-202, 1999. 

[6] R. J. McCann, "Polar factorization of maps on Riemannian manifolds", 
Geom. Funct. Anal, vol. 11, no. 3, pp. 589-608, 2001. 

[7] J.-D. Benamou and Y. Brenier, "A computational fluid mechanics solution 
to the Monge-Kantorovich mass transfer problem", Numer. Math., vol. 84, 
no. 3, pp. 375-393, 2000. 

[8] S. Angenent, S. Haker, and A. Tannenbaum, "Minimizing flows for the 
Monge-Kantorovich problem", SLAM J. Math. Anal, vol. 35, no. 1, pp. 61-97, 
2003. 

[9] G. Loeper and F. Rapetti, "Numerical solution of the Monge-Ampere equa- 
tion by a Newton's algorithm", C. R. Math. Acad. Sci. Paris, vol. 340, no. 4, 
pp. 319-324, 2005. 

[10] G. Carlier, a. Galichon, and F. Santambrogio, "From Knothe's transport 
to Brenier's map and a continuation method for optimal transport", SLAM J. 
Math. Anal, vol. 41, no. 6, pp. 2554-2576, 2009/10. 

[11] M. Rosenblatt, "Remarks on a multivariate transformation", Ann. Math. 
Statistics, vol. 23, pp. 470-472, 1952. 

[12] H. Knothe, "Contributions to the theory of convex bodies", Michigan Math. 
J., vol. 4, pp. 39-52, 1957. 

[13] A. FiGALLi, F. Maggi, and a. Pratelli, "a mass transportation approach 
to quantitative isoperimetric inequalities", Lnvent. Math., vol. 182, no. 1, 
pp. 167-211, 2010. 

[14] L. A. Caffarelli, "The regularity of mappings with a convex potential", J. 
Amen Math. Soc, vol. 5, no. 1, pp. 99-104, 1992. 

[15] A. FiGALLi, Y.-H. Kim, and R. J. McCann, Regularity of optimal transport 
maps on multiple products of spheres. To appear in J. Eur. Math. Soc. (JEMS), 
2011. 

[16] D. GiLBARG AND N. S. Trudinger, Elliptic partial differential equations of 
second order, vol. 224 of Grundlehren der Mathematischen Wissenschaften, 
Springer-Verlag, Berlin, 1977. 

[17] R. S. Hamilton, "The inverse function theorem of Nash and Moser", Bull 
Amen Math. Soc. (N.S.), vol. 7, no. 1, pp. 65-222, 1982. 



29 



